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ABSTRACT 

Context. Star-forming galaxies are a highly biased tracer of the underlying dark matter density field. Their clustering can be studied 
through the cosmic infrared background anisotropies. These anisotropies have been measured from 100 yum to 2 mm in the last few 
years. 

Aims. In this paper, we present a fully parametric model allowing a joint analysis of these recent observations. 

Methods. In order to develop a coherent model at various wavelengths, we rely on two building blocks. The first one is a parametric 

model that describes the redshift evolution of the luminosity function of star-forming galaxies. It was introduced in Bethermin et al. 

(2011) and compares favorably to measured differential number counts and luminosity functions. The second one is a halo model 

based description of the clustering of galaxies. Starting from a fiducial model, we investigate parameter degeneracies using a Fisher 

analysis. We then discuss how halo of different mass and redshift, how LIRGs and ULIRGs, contribute to the CIB angular power 

spectra. 

Results. From the Fisher analysis, we conclude that we cannot constrain the parameters of the model of evolution of galaxies using 
clustering data only. The use of combined data of Ce, counts and luminosity functions improves slightly the constraints but does not 
remove any degeneracies. On the contrary, the measurement of the anisotropies allows us to set interesting constraints on the halo 
model parameters, even if some strong degeneracies remain. Using our fiducial model, we establish that the 1-halo and 2-halo terms 
are not sensitive to the same mass regime. We also illustrate how the 1-halo term can be misinterpreted with the Poisson noise term. 
Conclusions. We present a new model of the clustering of infrared galaxies. Our framework allows a coherent and joint analysis of 
various probes of infrared galaxies : number counts, luminosity functions and clustering measurements. However such a model has a 
few limitations, as the parameters of the halo occupation suffer from strong degeneracies. 
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1. Introduction 



Infrared star-forming galaxies are mainly characterized by a very 
high star formation rate, tens or even hundreds times larger than 
that of the Milky Way, -10-100 M Q /year ( |Kennicutt|19"98] l. The 
large number of young stars are embedded in dust that absorbs 
UV starlight and reemits it in the infrared (IR), from 5 fj.m to 
1 mm. As a result, infrared star-forming galaxies emit most of 
their energy (-95%) in the IR. In the far IR (A > 200/itn) and 
sub-millimeter, observations are limited by confusion, small spa- 
tial scales are lost because of the poor angular resolution of the 
instruments. Thus we observe the Cosmic Infrared Background 
(CIB) ( |Puget et al.|1996[|Ftxsen et al.|19 98) which is the contri- 
bution of infrared star-forming galaxies integrated over the age 
of the Universe, and its anisotropies. In the near and mid IR a 
large fraction of the CIB is resolved into sources whereas in 
the far IR only a few percents is. These fractions can be im- 
proved thanks to the use of statistical methods. For instance, at 
long wavelength, |Qliver et al. ( 20T0| > directly resolved 15%, 10% 
and 6% of th e CIB at 250, 350 and 500 jum in Herschel/SPIRE 
data wh ereas |Papovich et al] ( |2004| l resolved 70% of the CIB at 
24 yum. |Berta et al.|(|2010|> resolved 45% and 52% of the CIB 



these fractions to 50% and 75% respectively. As a result, sources 
detected at 24 /mi are the main sources of the CIB around its 
peak at ~200 /urn. The CIB is dominated by objects that get 
more and more massive as the redshift increases from lumi- 
nous IR galaxies (LIRGs) with 10 n L Q < L !R < 1O 12 L at 
0.8 < z < 1.2 with intermediate mass, to extreme LIRGs with 
10 12 L o < L f fl < 10 14 £fr> that d ominate at z > 2 and with masses 



> 10 n M Q (|Caputi et al.|2006| . 



Magliocchetti et al. (20071 and Farrah et al. (2006i derived the 



at 100 and 160 jum respectively by integrating number counts 
in Herschel/PACS data. Stacking 24 /jm sources they increased 



two-pomt correlation function of Ultra LIRGs atz—i .6-2.7 and 
1 .5 < z < 3 respectively. They brought to light the very strong 
clustering of infrared star-forming galaxi es and their embed- 
ment in very massive halos of 1O I3 M . Cooray et al. ( 2010} 
computed the angular correlation function with Herschel/SPIRE 
data. They found that 250 /vm sources are in DM halos with 
masses around 10 12 M o that lie at z ~ 2.1 whereas bright 500 //m 
sources are in more massive halos ^ 10 I3 M o at z ~ 2.6. More 
recently, Magliocchetti et al. (2011 1 derived the 3D correla- 
tion function of infrared sources using Herschel/PACS data up 
to a redshift of 3. They obtain that their galaxies lie in haloes 
with > 1O 124 M , value that is in agreement with previous stud- 
ies. However the two-point correlation function is not easily 
computed using IR data because of confusion. As said before, 
confusion can be circumvented through the use of statistical 
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methods. Indeed, clustering can be measured in the correlated 
CIB anisotropies (CIBA). It has first been detected as an ex- 
cess of signal at intermediate scales by |Lagache et al. ( 2007 ) 
and |Grossan & Smoot| ( p007] l at 160 /mi in the Spitzer Multi- 
band Imaging Photometer (MIPS) data. These measurements 
have been followed by the detection in the Balloon-borne Large 
Aperture Sub-millimeter Telescope (BLAST) data at 250, 350 
and 500 //m ( |Viero et aL|2009|) a nd by that of the South Pole 
Telesc ope team (|Hall et aL 2010 1 at 1.3 and 2 mm. More re- 
cently Penin et al. ( 201 l| l measured the clustering signal by re- 
moving accurately the cirrus contamination at 100 and 160 /mi. 
The power spectrum of the CIBA has also been computed using 
Herschel/SPIRE at 250, 350 and 500 /mi ( jAmblard et al.|2011) , 
taking advantage of its angular resolution, and using Planck/HFI 
at 350 , 550, 850 /mi and 1.3 mm ( |Planck Collaboration et al.| 
201 1 J > taking advantage of its sky coverage. Therefore, the clus- 



tering of infrared star-forming galaxies in the CIBA has been de- 
tected over a large range of wavelengths and angular scales. All 
these results have been analyzed in several ways, hardly com- 
parable. As a first analysis, Lagache et al. ( 2007 1 derived the 
linear bias, the proportionality coefficient between the fluctua- 
tions of the dark matter (DM) density field and emissivities of 
galaxies. They found b = 2.4 + 0.2 and Viero et al. (2009) found 
b — 3 + 0.3 which implies that these galaxies are a highly biased 
tracer of DM. The difference between these two biases may be 
explained by the fact that at longer wavelength, higher redshift 
infrared star-forming galaxies are probed (Lagache et al. 2005; 
|Fernandez-Conde et al.[2 008) and thus are found to be more bi 



ased. New measurements needed more complex models. Viero 
|Amblard etaE] ( [201 1| > and |Planck Collaboration 



et al. 


(2009 


et al. 


(2011 



introduced a halo occupation distribution for the 
study of CIBA. It describes the DM distribution and especially 
how galaxies are distributed in one DM halo. It appears that each 
wavelength must be fitted separately which indicates an evo- 



lution of the clustering with the redshift ( |Planck Collaboration 
|et al.|20lT) . Most of the models determined the mass of the ha- 
los where infrared star-forming galaxies lie and thus where star 
formation occurs. 

In the long term purpose of analysing all these new measure- 
ments in a consistent way, we present a new model of the cluster- 
ing in CIBA. We use the halo model formalism ( |Cooray & Sheth 



(2002 )) which has been often used in the last few years to predict 



and to interpret galaxy clustering. We link it to a recent model 
of infrared star-forming galaxies evolution that reproduces well 
number counts and luminosity functions (Bether min et al.|201 1 



This model of clustering has been successfully used to fit Planck 
data ( |Planck Collaboration et al.|201 1) . 

The paper is organized as follows. We describe the model and 
its parameters in Sect. 2. We set a fiducial model inspired from 
Viero et al.|(|2009) ; |Planck Collaboration et al.| ( |20TT] l; |Amblard| 



et al. ( 201 1[ ) and compute angular power spectra for several in- 



struments with which we carry out a Fisher analysis in Sect. 4. 
Sect. 5 is dedicated to interpreting measurements such as the 
redshift and halo-masses contribution to the power spectrum, the 
linear bias, the influence of the mean emissivities and the con- 
tribution of LIRGs and ULIRGs to power spectra. We finally 
conclude in Sect. 6. Throughout this study we use the Wilkinson 
Microwave Anisotropy Probe 7-year Cosmology ( |Larson et al.| 
[20TT) . 

2. Why a new model? 

As said previously several models of clustering in the Cosmic 
Infrared Background already exist so why constructing a new 



one? There have been several measurements of the clustering 
and different models have been applied to analyze each mea- 
surement. Moreover their approaches are different which make 
comparison of the results difficult if not impossible. Thus one 
single model that ties together all available measurements is ap- 
pealing, especially to analyze them simultaneously and search 
for an evolution of the clustering. Such a model requires three 
ingredients: a DM distribution, a relation between galaxies and 
DM halos and an evolution of infrared star- forming galaxies. 



Lagache et aL] ( |2007| l, |Amblard & Coor ay (2007) and| Viero etal. 
(2009) used the model of galaxy evolution of |Lagache et al 



( |2003| l. This model was the most up to date model at that time. 
It fitted well differential number counts and luminosity func- 
tions measurements from 24 to 850 fim (IRAS, Spitzer/SCUBA). 
However it does not reproduce very well new measurements, es- 
pecially differential number counts from Herschel. This is a phe- 
nomenological model in which the evolution of the luminosity 
function was tuned to reproduce the constraints available at that 
time. It over-predicts the luminosity density at high-z. Moreover, 
it does not reproduce very well the observed redshift distribution 
of the CIB ( Jauzac et al.|201 l) . It predicts a peak at z ~ 1 that 
is not observed. The angular power spectra of CIBA strongly 
depend on the redshift distribution of the sources through the 
emissivities (see Sect. |3.2) , Therefore, a 'valid' distribution in 
redshift is important and a more robust model in agreement with 
most recent measurements is needed. 

|Hall et aL] (2010i used the galaxy templates from [Lagache et aL] 
( |2003[ l in order to check a simple model with a 'single SED'. 
This model has only a few parameters that can be changed eas- 
ily and thus adapted to each of their wavelengths. They fixed the 
shape of the power spectrum and only changed its amplitude de- 
pending on the wavelength. Amblard & Cooray (2007 ) used the 
same model of infrared galaxies evolution : they used the lumi- 
nosity function as a function of redshift at 350 /mi coming from 



Lagache et al. ( 2003| l that they matched to conditional luminos- 



ity functions (CLFs). Other wavelengths are extrapolated from 
the 350 /mi. Finally Amblard et al. 201 1 avoid using any model 
of galaxies evolution by letting free the redshift distribution of 
the cumulative flux coming from the background faint galaxies 
in several redshift bins. 

Concerning the distribution of DM, Viero et al. ( 2009 ) and 
Amb lard et al.| ( |20TT| ) used the formalism of the halo model and 
the same halo occupation number whereas Lagache et al. (2007) 
and |Hall et aT7| (|20 1 0)) considere d a linear power spectrum for 
dark matter. Amblard & Cooray (2007 ) also used the halo model 
formalism through CLFs. By integrating CLFs on the luminos- 
ity, the halo mass function is recovered. However this approach 
depends on too many parameters that cannot be constrained si- 
multaneously. 

In order to construct a new model, we link an up to date model of 
galaxies evolution to a recent version of the halo model. We use 
the model of evolution of galaxies of Bethermin et al. ( 2011) . 
It reproduces well Herschel measurements as well as older ones 
(from 15 //m to 1.1 mm). It also very well reproduces the red- 
shift distribution of the CIB of Jauzac et al. ( 201 1 1. We use an 
updated version of the halo model of Viero et aTji 2009), the halo 



occupation distribution (HOD) introduced by |Tinker & Wetzel| 
(2010). This HOD reproduces well the angular correlation func- 
tion of optical galaxies, red (star-forming) and blue (quiescent) 
galaxies at 0.4 < z < 2. Therefore we make a strong assumption 
here, assuming that this description would work on star-forming 
galaxies. Given the current lack of understanding of the details of 
the process of star-formation and its evolution with redshift, it is 
difficult to define what would be a better HOD prescription and 
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we therefore stay with this one. We study power spectra com- 
ing from our model for several wavelengths/instruments: 100 pm 
IRAS, 160 pm Spitzer/MIPS, 250, 350, 500 pm Herschel/SPIRE 
and 850 pm, 1.3 and 2mm Planck/HFI. A list of the available 
data of CIBA power spectra is given in table [2] 

3. The model 

3.1. The parametric model of star-forming galaxies evolution 

To reproduce the angular power spectrum of the CIBA we need a 
model for the redshift evolution of star-forming galaxies. We use 
the model presented in Bethermin et al. (201 1 1. It is a backward 
evolution model based on a parametrized luminosity function 
and on galaxies spectral energy distribution templates. 
Beth ermin et al.| ( [201 1) consider a luminosity function (LF) that 
behaves l ike a power law for L « L* and like a Gaussian for 
L » L* ( |Saunders etal.|1990| l : 



<5(L /fi ) 



dN(L IR ) 
dVd\og m (L IR ) 

. l-cr 



= $*(z) 



exp 



'2^ l0g ' 



1 + 



L*(z) 



(1) 
(2) 



where <t>(X/fi) is the number of galaxies with the infrared bolo- 
metric luminosity L; R within the comoving volume dV and the 
bin d\og m L. <£* is a normalization constant that fixes the den- 
sity of sources. The low and high luminosity parts have different 
slopes, I - a and I - a - 1/cr 2 / ln 2 (10) respectively. L* repre- 
sents the luminosity at the break. The parameters that describe 
the luminosity function are listed in table [T] 
The luminosity function evolves with the redshift through L* 
and <p* : 



Viz) = L*(z = 0)(l+z) ri 
0*(z) = </>*(z = 0)(l+z) r * 



(3) 
(4) 



Exponents ri and are not identical for all z. Two breaks are 
imposed to reproduce the evolution of the LF. The first one Zbreak 
is a free parameter and is found to be around 1 . The second one 
is fixed at z — 2 to avoid divergence at high z. Between these 
two breaks, the values of r L and change as shown on table [T| 



Bethermin et al. ( 201 1} used the SED library of |L"agache 



et al.| (|2004[). It contains two galaxy populations: star-forming 



and late-type galaxies. The latter emit half or less of their en- 
ergy in the IR whereas the former emit more than 95 % of their 
energy in the IR. The fraction of each population depends on lu- 
minosity. Indeed, late-type dominate at low luminosity whereas 
star-forming dominate at high luminosity. For a given bolometric 
luminosity, the fraction of star-forming is : 



fsF - 



®SF 

<D 

1 + tanh[o- pop (L)log n) (L IR /L pop )] 



(5) 
(6) 



L pop is the luminosity where <S>sf — ^late-npe and cr pop charac- 
terizes the width of the transition between the two populations. 
Differential number counts are then derived for each population 
and then summed. At flux S : 



dN . 
dS 



(S) 



a 



dN„ 



dlog w L IR dV 



pop 



dlog lQ L IR dV dS dz 



dz 



(7) 



02 




Fig. 1. CIB per redshift bins from the model of |Bethermin et al. 



(J20TTJ. The high redshift contribution increases with the wave- 
length. 



where dN/dS is the number of sources per flux unit in a unit solid 
angle and pop = late-type or pop — star-forming. 
The model of Bethermin et al. ( 201 \\ is described by thirteen 
free parameters. Best fit parameters and confidence areas are 
computed using Monte-Carlo Markov Chains on available and 
trustworthy differential number counts and luminosity functions 
at certain wavelengths. Bethermin et al. ( |2011[) fitted number 



count s of Spitzer/MIPS at 24, 70 and 160 pm ( |Bethermin et al 
2010b , those of Herschel/SPIRE at 250, 350 and 500 pm ([Oliver 
2010 ) and those of Aztec at 1.1 mm ( |Austermann et al 



et al. 



2010[ Scott et al.|2010 1. A couple of luminosity functions at dif- 
ferent redshifts are also fitted, the 8 pm one at z = 2 from[Caputi 
et al. (2007 1, others derived from Rodighiero et al. 2009 (a lo- 
cal LF at 24 pm, a 15 pm one at z = 0.6 and a 12 pm LF at 
Z = 1) and that at 60 pm at z = from Saund erset al.| l |1990) . 
Moreover, absolute measurements of the CIB are also used as a 
model constraint (Laga che et al.|19 99 ). We do not describe and 
discuss all the fits here (for a full discussion see |Bethermin et aL] 
|201 1[ ). Using the best fit, this model also provides the redshift 
distribution of the CIB as shown on Fig.[T We see that higher- 
redshifts contribution increases with wavelength : the redshift 
slice 0.3 < z < 1 dominates up to 400 pm whereas in the sub- 
millimeter z > 2 dominates. This model provides a very good 
agreement with the CIB redshift distribution (|Jauzac et al. 201 1 



We will study how these thirteen free parameters can be con- 
strained with power spectra of the CIBA. 

3.2. The angular power spectrum 

According to Haiman & Knox ( 2000} , Knox et al. (2001 1 and 
using the Limber approximation, the angular power spectrum of 
the anisotropies of the CIB at wavelengths A and A! is : 



dz 



dr a (z) 
dz~dT 



Mz)jAz)P ! 



s (k=-,z) 
d A 



(8) 



where I is the multipole, r is the conformal distance from the 
observer, a(z) the scale factor, d A the comoving angular diam- 
eter distance, and jx (z) the mean emissivity per comoving unit 
volume at wavelength A as a function of z. When A — A' we 
recover the auto power spectrum. P ss (k) is the galaxy three di- 
mensional power spectrum. 
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The emissivities are computed using the parametric luminosity 
functions following : 



where dN/dM is the halo mass function. We will use the uni- 



jv(z) 



I — 

dz 



fl 



S(L IR ,z) 



dN 



dzd(lnL IR ) 



d(lnLin) 



(9) 



where dN J 'dzd(lnL) is the number of galaxies per redshift bin 
dz and per luminosity bin d(lnL) and S the flux. Each galaxy 
population (late-type and star-forming) emissivity are computed 
and summed to get the overall emissivity. Fig.[2]shows emissiv- 
ities as a function of redshift. The two discontinuities at z ~ 0.9 
and z — 2 are due to the breaks imposed by the parametrization 
of the model of galaxies. It is clear that as the wavelength in- 
creases, the contribution from the high redshift part increases. 
Emissivities are color corrected according to their instrument 
and wavelengths to give Q in Jy 2 /sr (for the photometric con- 
vention v/ v =cst). 

In the context of the halo model, P ss (k) is the sum of the clus- 
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Fig. 2. Emissivities versus the redshift for different wavelengths. 
The contribution from high redshift increases with the wave- 
length. 



tering in one single halo (lh) and in two different halos (2h): 



Pss(k) = Pi h (k) + P 2h (k) 



where 



Pih(k) 



- f 

Jm 



dN {N ga i{N ga i- 1)> 



Im dM 
Pihik) = P lin {k) 



n i 

gal 

f —b(M) 
Jm dM 



U(k, M) p dM 

(Ngai) , 



n ga i 



-U(k, M)dM 



(10) 

(11) 

(12) 
(13) 



Here M is the halo mass, Pu„(k) is the dark matter linear power 
spectrum (computed with the fit of Eisenstein & Hu 1998), 
U (k, M) the normalized Fourier transform of the halo density 
profile that is assumed to be that of |Navarro et al. (T996) trun- 
cated at the virial radius. b(M) is the halo bias, {N ga i) the prob- 
ability of having N ga i galaxies in a halo of mass M and we con- 
sider p — 2 ( |Cooray &~S heth 2002). The mean number density 
of galaxies n ga i is given by: 



Igal 



versal form given by Tinker et al. ( 2008 ) as well as its redshift 
evolution. We use its associated halo bias (see Eq. Al in |Tinker| 
|eTaT1 ( |2009l )). 

The halo occupation number introduces galaxies in the halos sta- 
tistically. Recent data and simulations suggest a necessary dis- 
tinction between the major galaxy that lies at the center of the 
halo and the satellite galaxies that populate the rest of the halo. 
Above a given mass threshold, most halos will host a central 
galaxy. Above a second higher mass threshold, they will also 
host satellite galaxies. N ga i can thus be written as : 



(Ngai) = (Nen) + (N sat ) 



(15) 



According to the prescription of Tinker & Wetzel (2010 1, the 
occupation function of central galaxies is: 



_ 1 



1 +erf 



logM-logM m 



(16) 



where M m ,„ is the halo mass at which a halo has a 50 % prob- 
ability of hosting a central galaxy. cr log M controls the width of 
the transition between zero and one central galaxy. There is a 
smooth transition between low mass halos that do not contain 
bright enough galaxies to be seen in the data (M << M min ) and 
more massive ones that always contain a bright central galaxy. 
(M » M,„i„). The satellite occupation function is: 



(N sal ) = 



1 



1 +erf 



log M - log 2M„ 

V \os>M 



M 

M sm 



(17) 



It has a cut-off of the same form as the central occupation with 
a transition mass twice larger than that of the central to prevent 
halos which have a low probability of hosting a central galaxy to 
contain satellite galaxies. The number of satellite galaxies grows 
with a slope of a sat . Both number of galaxies as well as their sum 
are plotted on Fig. [3] 

With this model, the angular power spectrum of CIBA depends 
on only four halo model parameters a sat , M m j„, M sa , and cr/ ogM . 
Cosmology is fixed at WMAP7 values. Our parameters are listed 
in Table [T] with their meaning and their fiducial values that we 
set in Sectg]. 

The long term purpose of our model is to look for best fits of 
these parameters for Spitzer/MIPS, IRIS, Planck, Herschel and 
SPT data and study their evolution with wavelength. However it 
is beyond the scope of this paper. Therefore we will not com- 
pare the data to the power spectra coming from our model. 
Our first aim here is to study the parameter space and to inves- 
tigate particularly the behaviour of the halo bias, the halo mass- 
contribution to the power spectrum and its redshift distribution. 
To do so we will consider a set of fiducial halo parameters iden- 
tical at all wavelengths. 

4. Power spectra and parameters degeneracies 

In this section, we present the CIB power spectra computed 
with the model detailed in the previous section for several wave- 
lengths in the far-IR and sumillimeter. We then study the de- 
generacies of the parameters, looking first at the galaxies model 
parameters and second at the HOD parameters. 



4.1. Power spectra 

Our fiducial model is 
logM„„„ = 11.5, M mt = 



set with the HOD parameters 
10M,„,„ and a = 1.4 at all wave- 



(1 ) lengths. These values are motivated by the parameters fit of 
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Fig. 4. CIB anisotropics power spectra at several wavelengths. The continous black line shows the power spectra of the clustering. 
The dotted black line is the 14ialo term of the power spectrum and the dashed line is the 24ialo term. The blue horizontal line 
represents the shot noise level and the red dot-dashed curve is the total power spectrum. Titles of the plots give the reference for the 
shot noise level and the used bandpass filters. 
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parameter name 


Description 


Value 




Minimal mass of a halo to have a central galaxy 


1U Mq 


M sa , 


Nomalisation mass for satellite galaxies 


inl2.5 tla 
1U Mq 




Slope of the number of satellite galaxies at high mass 


1 A 


0~ logM 


Scatter in halo mass 


n ia q 


a 


Faint end slope of the IR bolometric LF 


1 TT2 

l.zzj 


(T 


Parameter driving the bright end slope 


U.4UO 


L*(z=0) 


Local characteristic luminosity of the LF 


2.377x10 L Q 


«f>*(z=0) 


Local characteristic density ol the LF 


3.234x10 gal/dex/Mpc 


r UJz 


Evolution ol the characteristic luminosity between and Zbreak 


2.931 


r <t>*Jz 


Evolution ol the characteristic density between and Zbreak 


0.774 


Zbreak 


Redshift of the first break 


U.o ly 


r Li,,mz 


Evolution of the characteristic luminosity between Zbreak <md 2 


A 737 
t +. I j 1 


r <!>-*, mz 


Evolution of the characteristic density between Zbreak an d 2 


-6.246 


r U,hz 


Evolution of the characteristic luminosity for z > 2 


0.145 




Evolution of the characteristic density for z > 2 


-0.919 


^pop 


Luminosity of the transition between normal and star-formig templates 


23.677x10'° 


& pop 


Width of the transition between normal and star-forming templates 


0.57 



Table 1. Parameters of our model. The first part of the table lists the halo model parameters and the second part lists the parameters 
of the model of galaxies. The values of the latter are the mean ones of Bethermin et al. (2011|l. 
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Fig. 3. Number of galaxies versus dark matter halo mass. The 
blue dashed line shows the central galaxies, the red dotted- 
dashed line shows satellite galaxies and the black continuous 
line shows the total. We use the parameters of our fiducial model 
(see Sect. |4|>, that is to say logM,„,„ = 11.5, M sat = 10M„„„ and 
a =1.4 



Viero et al.|(l200^ , [Pranck Collaboration et al.| ( |20TT) , |Amblard| 
et al.| (201 1 1. As the halo parameters slightly depend on the 
wavelength (in reality but it is not the case here), the power 
spectra presented in this section may not be seen as an exact 
prediction but as a basis for a qualitative study. For this fiducial 
model, we present on Fig. [4] the power spectra for different 
experiments and selected wavelengths, from 100 fim to 2 mm. 

The comparison to measurements also requires the introduc- 
tion of a shot noise term due to the finite number of galaxies. We 
compute it using our galaxy evolution model (Bethermin et al. 
[20TT) . 



Lcl us 



+ c 



[.shot 



(18) 



where Q jC ; M is the power spectrum of the clustering and C( tS h t 
is the shot noise. C^shot depends on the flux cut applied to the 



data when removing or masking the brightest sources. Typical 
flux cuts for different far-IR and sub-millimeter experiments are 
given in Tab. [2] On Fig. [4] we only show one shot noise level 
per wavelength corresponding to the instrument given in the plot 
titles, for purpose of clarity. 

The contribution of the 1-halo (2-halo) corresponds to the dotted 
line (dashed line). The instrument noise is not shown here but it 
is often negligible, the CIB being measured with a very high S/N 
even at spatial scales close to the angular resolution. Looking at 
this plot, we see clearly that the interplay between 2h, lh and 
shot noise terms will make the interpretation of measurements 
quite subtle. The contribution of the 1-halo term decreases with 
the wavelength which can prevent its measurement if the reso- 
lution of the instrument is low. For example, Planck with its 5 
'angular resolution at high frequency, cannot measure directly 
the shot noise level and the 1-halo term can be easily misinter- 
preted as shot noise. Reversely, the 1-halo term dominates a large 
range of scales at 100 and 160 /mi and thus can be measured ac- 
curately at these wavelengths. 

High wavelengths probe higher redshifts than short ones and ha- 
los are relatively smaller at high redshifts compared to those in 
the local Universe as the latter had time to accrete more matter. 
Therefore the scale of the intersection between the 1- and 2-halo 
terms shifts towards higher ( as the wavelength increases. It goes 
from I ~ 50 at 100 fimtoi ~ 1000 at 2 mm. | Viero etal] (2009) 
also observed this trend. At 250 yum the crossing is at k ~ 0.03 
arcmin -1 {I ~ 648) whereas it is at k ~ 0.06 arcmin -1 (£ ~ 1296) 
at 500 jjm. The exact crossing point differs from ours because of 
the HOD parameterization. 



4.2. Variation of power spectra with the galaxy-evolution 
model parameters 

In order to do an overall study of how our model parameters 
can be constrained, i.e. to investigate how degenerated they are, 
we construct the Fisher matrix associated to the power spectra. 
We write the Fisher matrix for angular power spectrum measure- 
ments as: 



1 dC x ( dC\ 



(19) 
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wavelength (pm) Instrument 



Reference 



Flux Cut (mjy) Shot noise level (Jy 2 /sr) 



100 



160 

^50~ 



350 
350 



500 
550 



850 



1363 
1363 
1380 



2000 
2000 
2097 



IRIS 



Penin etal. J20111 



Spitzer/MIPS 
Herschel/SPIRE 



Lagache et al. ( 2007 ) 
~Amblardetal.~(201ir 



Herschel/SPIRE 
Planck/HFI 



Amblardetal. (20111 



Herschel/SPIRE 
Planck/HFI 



Planck Collaboration et al 



[Amblard et al.||20TI 
Planck Collaboration et al. ( 201 1 1 



Planck/HFI 



Planck Collaboration et al 



SPT 
ACT 
Planck/HFI 



SPT 
ACT 
Planck/HFI 



Hall et a l. (2010) 
Fowler & Atacama Cosmology telesc ope Team 1 2010| 
'Planck Collaboration et al. | poll) 



|Hall et al.||2010[> 

Fowler & Atacama Cosmology Telescope Team 1 2010| 



700 



200 

^0~ 



50 
710 



50 
540 



325 



6.4 
20 
160 



6.4 
20 
245 



7364±1232 



10834±3124 
6715+1458 



4362±1250 
5923±367 



1156±434 
1150+92 



138+22 



11.9±4.0 
12.5+3.9 
12.9±2.9 



1.73±0.54 
1.78±0.60 
1.4±0.3 



Table 2. Shot noise levels in Jy /sr from Bethermin et al. (2011 1 for available measurements of CIBA power spectra 



wavelength (/im) 


Instrument 


Jy7srto//K^. MB 


100 


IRAS 


9.59xl0 2i 


160 


Spitzer/MIPS 


3.12x10" 


250 


Herschel/SPIRE 


1.34 xlO 3 


350 


Herschel/SPIRE 


2.78 xlO" 1 


350 


Planck/HFI 


2.00 xlO" 1 


500 


Herschel/SPIRE 


7.45 xl0~ 4 


550 


Planck/HFI 


2.94 xl0~ 4 


850 


Planck/HFI 


1.20 xlO" 5 


1380 


SPT 


4.39 XlO -6 


1380 


Planck/HFI 


4.32 xl0~ 6 


2000 


SPT 


6.10 xlO -6 


2097 


Planck/HFI 


7.31 xl0~ 6 



Table 3. Conversion factors from Jy 2 /sr to /uK CMB . One should 
multiply the power spectrum in Jy 2 /sr (with the convention 
vl y =cst) by the coefficient to get jjK^ MB . 



where 07 are the errors on the measurements and they include 
both the cosmic variance and the instrumental noise at a multi- 
pole I: 



a-, = 



c + N <* 

W H t 

R 1 
D t ) 



fskyW + l) 



(20) 



where is the fraction of the sky we consider, N[ is the 
level of the instrumental noise and B 2 the power spectrum of 
the beam. To compute the Fisher matrices, we generate mock 
power spectra using our fiducial model and error bars derived 
following Eq. ( pO} , from 100 /urn to 1.3 mm. The range of multi- 
poles is taken to be consistent with the available data. At 350 
and 550 fim, we assume combined power spectra for Planck 
and Herschel, and we thus extend Planck power spectra to the 
Herschel limit in multipole. We plot Q in Jy 2 /sr. They can be 
converted in yuK 2 using the coefficients given in Tab. [3] 
The bottom left panel of Fig.|5]shows confidence ellipses (lcr in 
green and 2<x in black) coming from Q when trying to measure 
only the galaxy model parameters ri*^z an d r L*,iz- Clearly, they 
are very poorly constrained. For instance ri*j lz = 0.145± 15.55, 
or r L*,iz - 2.93±20.0. For reference, the constraints obtained us- 
ing current number counts are ri*^ z = 0.145+1.05 and ri*j z = 
2.93±0.27, as shown on the top left panel of Fig. [5] The fact is 
that these parameters enter in the expression of the Q through 
the emissivities which are integrated on all redshifts, they are 
thus hard to measure from clustering measurements alone. 



As a matter of fact, the lack of information in Q partially 
comes from the large number of parameters in the model of evo- 
lution of galaxies. To quickly quantify this we vary only a few 
of these parameters {r Li , Jz , r Li , M , z break , r M and r^ M ) assum- 
ing that the others are perfectly known. Fixing all but these pa- 
rameters corresponds to assuming that only the redshift evolu- 
tion of the LF is unknown, clearly an irrealistic assumption. Not 
surprisingly, while some of the degeneracies remain in these re- 
duced parameter space, on the whole, parameters are better con- 
trained. For instance, we now obtain ri+j z = 2.93+0.10 and rj*^ 
= 0.145+7.05 which are about a two order of magnitude and a 
factor 2 improvement, respectively as compared to the numbers 
above. 

In order to illustrate further this lack of information in Qs, we 
show how they change with only one parameter, ri*j z . We make 



it vary by +2cr from its best fit (<r coming from Bethermin et al. 
201 \\ . This parameter governs the evolution of the luminosity 
function for < z < Zbreak- A higher ri*,i z means a faster in- 
crease of the luminosity, thus a higher value of L*(z = Zbreak)- 
The top panel of Fig.[6]shows the influence of this parameter on 
the counts at 160 fim. A higher (smaller respectively) r L * j z im- 
plies higher (smaller) number counts thus more (less) galaxies 
on a large range of fluxes. This leads to a higher (smaller) emis- 
sivities as shows in the second and third panels of Fig. [6] This 
results in a modification of ~ 20% on the emissivities and from 
15 to 35 % on power spectra depending on the scale. The fact 
that this ratio is not constant is due to to the fact that the ratio of 
the emissivities is not constant with redshift (3rd panel of Fig.|6]l. 
We can see that all power spectra are consistent within error bars 
and thus we can hardly discriminate between them. Therefore, it 
is hard to constrain the evolution model of galaxies using only 
power spectra. 

More relevant data are required. We compare our confidence el- 
lipses with those obtained with luminosity functions and number 
counts data. To do so we use the covariance matrix of Bethermin 
|et al.| ( |201 1[ ). The error bars are in general much smaller and 
there are only a few degeneracies. For instance, as shown on the 
left panel of Fig. [5] ri+j K and r^*^ are still strongly degenerate 
but they are now much better constrained. 
So far, we have investigated how galaxy evolution parameters are 
degenerated and constrained using LF/counts and C( separately. 
The next step is to look at the degeneracies when combining 
all these data. To do so, we add the two Fisher matrices com- 
ing from the counts/LF and the C(. The bottom panels of Fig. [5] 
show the confidence ellipses for r L * j lz , r L * j z and r^* using the 
combined data. Note that the axis scales are different. The con- 
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-3.540-1.896 0.058 1.702 2.3 2.7 3.2 3.6 
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Fig. 5. lcr (dashed green), 2cr (dash-dot blue), 3cr (black) likelihood contours of the galaxy evolution model parameters ri*j lz , 
rv,h and r^j K . Top, Middle and Bottom panels show the contours computed using counts/LF data, C( data only, and combined 
counts/LF and Q data, respectively. The continous line shows the direction of degeneracy using only counts/LF, the dotted line 
shows that using Q and the dashed line is that using all data. Q alone cannot constrain the galaxy-evolution model parameters and 
the combination of C( and counts/LF data do not change much the constraints. It slightly improves them. 



tinous/dashed/dotted lines on Fig. [5] indie ate the degeneracy di- 
rections. They are different which illustrate the complementarity 
of the two data-sets and the constraints can be greatly improved. 
For example, the errors on r^* /,, are decreased by a factor of 1.5 
but the errors on r L * j lz is not changing. However, this plot also 
clearly shows that overall, the number counts and LF measure- 
ments are much more powerful when looking at constraining the 
LF. However, Q can still constrain the global evolution of galax- 



ies through their mean emissivities. A first attempt was made by 
Amblard et al. (2011 ) who did not use a model of galaxies to 



compute the emissivities but bin them in several redshift inter- 
vals and considered the values of the emissivities in these four 
bins as free parameters. They also required that the integrated 
source density is within the 68% confidence level ranges of the 
CIB obtained by FIR AS. 
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multipole L 



multipole I 



multipole I 



Fig. 7. CIBA power spectrum at 160 fim obtained using several values of a sau M m i n and M saT . The level of the Poisson noise has also 
been added to the power spectra (black continous line). Pink dots are the data from Lagache et al. ( |2007 1 at 160 fim. When fixed, 
the parameters are those of the fiducial model, a sat — 1.4, M min = 10 11 5 M and M sat = W AJ M Q . Left panel: The blue line is for 
a sal = 1 .8, the green one for a sat = 1 .4 and the black one for a sat = 1 . Middle panel: The blue line is the clustering power spectrum 
for Mmin = 10 I2 5 M Q , the green one for M min = 10 1L5 M G and the black one for M min = 10 I0 5 M o . Right panel: The blue line is for 
M sa , = 10 13 - 5 M G , the green one for M sat = 10 12 5 M Q and the black one for M sat = 1O 1L5 M . As expected, Qs strongly depend on 
the halo parameters. We can expect strong degeneracies between those parameters. 



4.3. Halo occupation distribution parameters and their 
degeneracies 

The shape of the clustering power spectra strongly depends on 
the HOD parameters as shown on Fig. [7] We vary M mi „, M sal 
and a S at. Each panel shows the power spectrum at 160 fim mea- 
sured by Lagache et al. (2007) in addition to the ones coming 
from the model. In each plot we vary only one HOD parame- 
ter and hold the others fixed to the values of the fiducial model, 
a sat = 1.4, M mi „ = 10 1L5 M G and M m = 10 125 M G . Both the 
shape and amplitude strongly vary. The fact that some similar 
changes are observed using different parameters suggest strong 
degeneracies. We compute the Fisher matrix as in Sect. |4.2| and 
likelihood contours at 1- and 2-cr are shown on Fig. [8] The error 
bars on <t; os m are very large, thus we fix its value to 0.65, follow- 
ing |Tinker & Wetzel| ( |2010| l who studied the galaxy-clustering 
in optical surveys. Using Q only we are not able to constrain 
its value. M„„„ and M mt happen to be highly degenerated in the 
direction M sa , = 3.3M m j„. 

In previous works using optically selected galaxies, a is usually 
set to 1 ( Gao et al. ( 2004 1) and M„„„ and M sal are the only pa- 



rameters to be fitted to the data. Indeed, a is hardly constrained. 
For instance |Tinker & Wetzel| ( |2010[ ) used the same halo dis- 
tribution number on a sample of red and blue galaxies in the 
range 0.4 < z < 2. They fitted well their correlation function 
fixing a — 1 and letting free M mi „ and M sat . But when they let a 
free in addition to the two others, they obtain unrealistic values 
for a ( |Tinker & We tzel 2010). However, setting a — 1 might 
not be appropriate for CIBA. The halo model is commonly used 
in galaxies catalogs that are not deeper than z ~ 2 and CIBA 
probe higher redshifts, especially at long wavelengths. Moreover 
optically-selected galaxies are not only star-forming galaxies, 
and there is no reason why optically-selected galaxies and star- 
forming galaxies should behave in the same way. Magliocchetti 
et al. ( 2008[ ) used a similar form of the halo model to analyze the 
angular correlation function of 24 fim sources at 0.6 < z < 1.2 
and z > 1.6. Using two halo density profiles (NFW and a steeper 
one p ~ r -3 ), they derived a ~ 0.7 for the steeper profile and 
a ~ 0.8 for the NFW one. They get the same results for both 
sets of data. On the contrary, |Cooray et al.| ( |2010| computed the 
angular correlation function of sources detected at 250, 350 and 
500 fim in Herschel/SPIRE data. They used the same halo distri- 
bution as ours and get a - 1.3 ± 0.4, a < 1.8 and a < 1.6 at 250, 



350 and 500 ^m respectively. Finally, the Planck Collaboration 
|et al.| ( [20TT|) derived values of a compatible with 1. The discrep- 
ancy with lMagliocchettfe t al. (2008 ) may be due to the different 
properties of the bright galaxies selected at 24 fim and those that 
contribute to the CIBA at longer wavelengths. Here combining 
all CIBA measurements from 100 fim to 1.3 mm and low to high 
multipoles, a is well constrained and it is not strongly degener- 
ated with other parameters. Note that values of a > 1 implies 
that higher-mass halos contribute relatively more than smaller- 
mass ones compared to the halos in which lie optical galaxies 
such as those used by Tinker & Wetzel (2010). 
In our analysis of the degeneracies of the halo parameters we 
have only considered a set of parameters identical for all wave- 
lengths which is not the case in reality. It could thus be that the 
degeneracies depend on wavelength. We therefore checked that 
the degeneracies were not significantly changing when we com- 
pute the Fisher matrix with various set of parameters correspond- 
ing to the wavelength best-fit models. 

The halo parameters cannot be constrained by counts or LF as 
they only intervene in the clustering of galaxies in the equation 
of the C( (see Eq.|8]l. Therefore we cannot carry a joint analysis 
of the degeneracies of the halo parameters using counts/LF and 
C{ data all together. Note that in principle, we could extrapolate 
the number count measurements to constrain the total number of 
galaxies, which also depends on the HOD parameters, but this 
would be a difficult measurement as it would be strongly depen- 
dent on the flux cut, for example. 

Emissivities are given by the model of galaxy evolution but we 
want to investigate the degeneracies if they are binned in redshift 
and their values considered as free parameters as in |Am blard 
|et al.| ( p0TT) . They carried this analysis at redshift between and 
4. In order to be coherent with what has been done previously 
and to take advantage of our redshift range, from to 7, we split 
the whole redshift range in four bins, < z < 0.9, 0.9 < z < 2, 
2<z<3.5, 3.5<z<7. For each bin i we take the mean 
value of the emissivity that we call fry, with i = {1,2,3,4} and 
we compute the Q at 350 fim (we assume a combined Planck 
and Herschel power spectrum) and the associated Fisher ma- 
trix. Confidence levels are given on Fig. [9] First we see that 
the halo occupation number degeneracies do not change much 
(see the previous paragraph) apart from the error bars which are 
much larger. M mi „ is still strongly degenerated with M sat such 
as M sat - 2.1M min . The direction of the degeneracy is roughly 
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Fig. 9. lcr (green) and 2cr (black) likelihood contours of the halo model parameters computed with mock data at 350 jim. Instead 
of using emissivities coming from the model of galaxies, we split Q in four redshift bins, on which we use the mean value of 
the emissivity on each bin bp that we consider as free parameters. The bp are hardly constrained with our data. a sat is strongly 
degenerate with bj t \ and bjp which are the mean emissivities at low z. M m - m M mt and cri ogM are degenerate with bj$ and b u \, the 
mean emissivities at high z. 



the same as that derived using emissivities of the model as well 
as those of <Ti og M and M mt and M„„„. Therefore, the degeneracy 
directions are all similar, using the emissivities or letting them 
free. 

The degeneracies of the bj with the halo parameters depend on 
the redshift. Indeed, a sa , is highly degenerated with bj\ (0 < z < 
0.9) and b j2 (0.9 < z < 2) and not at all with bp (2 < z < 3.5) 
and bj4 (3.5 < z < 7). Therefore a sa , is constrained by z > 2 
galaxies whereas the other halo parameters behave in the oppo- 
site way, they are not degenerate with bj\ and bp but with bp 
and bj\. This redshift dependency is emphasized by their own 
degeneracies. bj\ and bp are strongly correlated, the same is true 
for bp and bu,. The degeneracies using other wavelengths are 
only slightly different, we do not show them here. In general, the 
couples (bp, bp) and (bp,bj4) are always strongly degenerate 
and the bp are degenerate with the halo parameters as shown on 
Fig. [9] Such a degeneracy dependence with the redshift has to be 
still understood. 



( 201 we carry a similar analysis using their redshift bins, that 
is to say, < z < 1, I < Z < 2, 2<z<3 and 3 < z < 4. We 
observe the same behavior as described previously but differ- 
ent from their results. They used Monte Carlo Markov Chains to 
compute the degeneracies and usually the two dimensional prob- 
ability distributions have two peaks (see their fig. S 13). Their S , 
are equivalent to our bp. S i is degenerate with the three oth- 
ers St, whereas there is no degeneracy between S 3 and S4. We 
do not discuss the degeneracy of the halo parameters with 5, 
as their parameterization of the halo occupation number slightly 
differ from ours. 

Using bp and fitting them on the data avoids us to rely on a 
model of evolution of galaxies. However they are poorly con- 
strained with the present data. Moreover the degeneracies be- 
tween the bp and the halo parameters strongly depend on the 
halo parameterization used. 



In order to compare our results with those of Amblard et al. 
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5. Interpreting measurements 

Now that we have determined a model and the associated param- 
eters degeneracies, we discuss their physical interpretation. 



5.1. Redshift and halo-masses contribution to the power 
spectrum 



The left panel of Fig. 15 shows the contribution to the Q by 
several redshift bins. As stated previously, the shorter the 
wavelength, the more important is the relative contribution from 
the low redshift. For example, while z < 0.7 contributes signif- 
icantly at 100 and 160 fim, it becomes much less important in 
the millimeter range. Reversely, the high redshift bin (z > 3) is 
negligible at short wavelength but has an increasing contribution 
when the wavelength increases. The redshift distribution can 
change with the choice of the halo parameters as shown on Fig. 
1 1 We have changed the values of a sat on the left panels and 
that of M„,i„ on the right panels (the change is in the lcr error 
bars of the best fit found by Planck Collaboration et al. (201 1 )), 
the other parameters are those of the fiducial model. We recover 
the trend noticed above more or less emphasized. 
Further, not only the redshift of the galaxies probed depends 
on the wavelength, but so does the mass of the halos in which 
they are embedded. Fig. 10 shows the contribution of mass 
and redshift to the 1- and 2-halo terms (at ( = 2002 and 
I = 100 respectively) from 100 fim to 2 mm. High-mass halos 
(M > 10 13 M Q ) contribute the more to the 1-halo term from 
100 fim to 2 mm at low redshift, and this dominant mass 
range stays constant with wavelength. More massive halos 
contain more galaxies than smaller ones therefore the galaxies 
contained in those halos contribute more to the angular power 
spectrum. This can be explained qualitatively in the following 
way. According to the mass function, at say z = 0.5, there are 
one thousand times more halos of M = l0 n M B /h than halos of 
mass M = lO 14 M //z. According to the HOD, one out of one 
hundred 10 n M o //i halo hosts a galaxy, whereas M - l0 14 M @ /h 
mass halos hosts on average 10 galaxies. Since the contribution 
to the 1-halo terms goes like N^ al , lore massive halos contribute 
relatively more to the 1-halo term. Note that this trend is less 
prononced for the 2-halo term since it goes like N ga /. At all 
wavelengths, as the redshift increases, the dominant mass range 
decreases to M ~ 10 n ~ 13 M Q as halos at higher redshifts are 
smaller than those at z — 0. 

The 2-halo term does not exhibit the same behavior : at short 
wavelengths, halos in a large range of mass 10 11 - 10 15 M o at 
low z contribute to the power spectrum. In parallel, intermediate 
masses contribute at higher z. As the wavelength increases, 
the relative contribution between high mass at low z and 
intermediate mass at high z becomes close to unity. It ends 
up in an equal contribution from high mass at low z and from 
intermediate mass at high z at 2 mm. Intermediate mass halos 
are more abundant which explains their high contribution. Both 
the 1 and 2-halo terms are sensitive to different mass regimes 
which evolve with the wavelength and thus with the redshift. 
Magliocchetti et al. (2008 ) selected 24 fim sources at 1.5 < z < 3 
and deduced from their correla tion function that they lie in 
10 I3 M o halos. |Vieroetal.| ( |2009ll found an M ef f ~ 10 1 3 - 2 M Q for 
unresolved galaxies at z > 1 at 250, 350 and 500 fim. Brodwin 



that LIRGs lie in halos M > 3 x 10 13 M G . All these results are 
in agreement with ours. Overall, halos with masses such as 
1O 12 ~ 13 M contribute the most to power spectra at all redshifts. 
However, we want to reemphasize here that these conclusions 
are model-dependent and depend on the particular emissivity 
model used, as discussed before. This is particularly true for the 
higher-z contribution, say z > 2. 



5.2. Linear bias 

Within our halo model we derive the linear bias as a function of 
the redshift following: 



bii„(k, z) 



(21) 



where P gg (k,z) is the galaxy-galaxy power spectrum coming 
from our model, Pu„{k,z) is the linear DM power spectrum and 
bii n {k,z) the linear bias. 

We will call the effective bias coming from the HOD model, 
i>HOD- On large scale uik — > 0, M) ~ 1, so the HOD bias from 
Eq.[l3]is: 



■>HOD 



(Z) 



I 



dN < Ngoi > 

dM b{M) s -^— 

dM n ga i 



(22) 



We plot in Fig. 12 the linear biases and the HOD bias for our 
fiducial model as a function of redshift for k = 1 (Mpc/h)' 1 
and k = 10~ 2 (Mpc//i) _1 where we also add current measure- 
ments detailed in Tab. [4] Note that in the linear regime, the HOD 
and linear biases are identical which is the case at k = 10~ 2 
(Mpc//z)~' but not at k = 1 (Mpc/Zz) -1 . On small spatial scales, 
the u(k — > 0, M) ~ 1 is not true therefore it is strongly different 
from the linear bias for k — 1 as shown on Fig. 12 



For both scales, our linear biases as well as the HOD bias is 
not in agreement with measurements. Neither of them show the 
same trend as the data points. The measured linear biases as well 
as the HOD biases grow quicker towards higher values than the 
biases extracted from our model. 



HOD bias measurements are from different HOD. Cooray et al. 
(2010 1 used the same halo occupation number to fit the corre- 
lation function and they found different parameters than ours. 
When using their parameters set, we do recover their results. 
The discrepancy could be explained by the fact that these mea- 
surements result from correlation function analysis, thus from 
resolved sources which is not the population we are studying 
here. 

Concerning the linear bias determined with unresolved galaxies 
( Lagache et al.|2007 Viero et al.|2009| l, our linear bias is not in 
agreement with the measurements either. Indeed when using un- 
resolved sources, the determination of the bias requires the use 
of emissivities, which are strongly model dependent (see Sect. 
5.3b and it can affect the bias. 



5.3. Influence of the mean emissivities 



et al. (20081 derived the angular autocorrelation function of (2004b. In Fig. 13 



dust-obscured galaxies selected with a color criterion. They 
determine d that they are in halos with an average mass of 



Previous models such as those of Laga che et al.| ( f2007|l and|Viero 
et al.| (|2009] > have used emissivities coming from |Lagache et aL 



we plo t the emissivity used in this paper as 
well as the Lagache et al. (2004) ones for reference. The peak at 
z ~ 1 in our emissivities is due to the parameterization of the LF. 



10 M Q . Gilli et al. (2007 1 selected star forming galaxies at 



24 //m, derived the projected correlation function and found 



Despite the shapes of the emissivities of Lagache et al. (2004) 
and ours are different, they display similar trends. The relative 
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wavelength 


kind of galaxies 




reference 


< z > 


bnoD 


bun 


24 


resolved 


Magliocchetti et al. (2008 1 


0.79 


1.70 




24 


resolved 


Mag 


liocchetti et al. (2008) 


2.02 


6.17 




24 


resolved 


Brodwin et al. ( 20U8 1 


2 




3.1-5.3 


24 


resolved 




Ciilh et al.|((10U7[> 


0.75 




1.3+0.1 


24 


resolved 




Uilh et al.|{z0D7} 


0.5 




1 A i Ci 1 

1.4±U.l 


9 /I 

24 


resolved 




Ulli et al.|jlDu7t 


1 

1 




l.o±U.Z 


100 


rpcnl vpH 

1 E.3UI V HI 


Saunders et al. ( 1V»92 ) 


o 




0.86 


160 


background 


Lagache et al. ( 2007 ) 


1 




2.4±0.2 


250-350-500 


background 




Viero et al.|^0U9| 


1.5 


2.2±0.2 


3+0.2 


250 


resolved 


c 


:ooray et al](|Z010|> 


9 i +0.4 
-0.7 


2.9±0.4 




350 


resolved 


|(Jooray et al.|<pfTU|> 


9 -3+0.4 


3.2±0.5 




250 


resolved 


|(Jooray et al.|(p010) 


26^ 

^•"-0.7 


3.6±0.8 





Table 4. Linear and effective bias measurements. The third column gives the mean redshift of the galaxies probed and the last but 
one lists the HOD/effective bias values and the last one gives the linear bias. 



contributions of high redshifts increases with wavelength while 
the contribution of low redshifts decreases. According to Jauzac 
et al. ( 2011| l, the model of Lagache et al. (2004) predicts too 
much power at high z. As the latter is forced to reproduce levels 
of the CIB and number counts, it does not predict enough power 
at low z. Therefore it predicts more galaxies at high z and less at 
low z. To illustrate how it influences our results, we show on Fig. 
14 the ratios of the power spectra computed with our emissivi- 
ties and those from Lagache et al. (2004). Up to 550 //m the ratio 
is around 1 up to I ~ 10000 and it increases strongly at higher I. 
At longer wavelength, the difference is much larger. It is in line 
with the over-prediction of power at high redshift of the model 
of Lagach e et al.| (|2004). The same halo parameters have been 
used for this plot, however it is clear that when fitting the model 
to the data with both emissivities we will not find the same halo 
parameters. 

As said previously the Lagache et al. (2004 ) model predicts too 
much power at high z, thus we need to compare the contribution 
in redshift to the Q. They are given on Fig. 15 In both cases 
we observe the same trend : high redshift contribute more and 



to their evolution with the wavelength. To do so we split the 
emissivities in the following way : 



more as the wavelength increases. However, with |Lagache et alT 
(2004) emissivities this evolution goes faster. For instance, at 
250 //m, the contribution of 1.5 < z < 3 galaxies is of the same 
order of magnitude that those of 0.7 < z < 1 .5 whereas in our 
case the former is more than one order of magnitude inferior to 
the latter. As we go towards longer wavelengths, the two highest 
redshift bins have an increasing contribution, and both dominate 
the power spectrum at 850 //m using Lagac he et al.| ([2"004 ) emis- 
sivities whereas using our emissivities 0.7 < z < 1.5 galaxies 
contribute also strongly. At 1.3 and 2 mm both the highest red- 
shift bins contribute the most but in our case, only the highest 
redshift bin dominate and the 1.5 < z < 3 bin has a smaller 
contribution. Therefore the shape of the emissivities strongly in- 
fluences our results, parameters determination and redshift dis- 
tribution. The interpretation of clustering measurements is thus 
based on the use of a reliable model of evolution of galaxies. 



5.4. Contribution of LIRGs and ULIRGs 

Star-forming galaxies are split in several categories according 
to their luminosities. Normal, luminous infrared, and ultra lu- 
minous infrared galaxies have luminosities L/« < 10 n M o , 
10 n M G < L IR < 10 I2 M Q , L m > 10 12 M G , respectively. LIRGs 
dominate the infrared energy output at z ~ 1 and ULIRGs at 
z ~ 2 (|Perez-Gonzalez et al.^OOSHLTFloc'h et al.|2005||Caputi 
|et al.| 2007 ) therefore we look at their contribution to the Q and 



; --normal . -;L1RC , ~jULlRG 
U = J A + J A + J A 



(23) 



and this contribution is squared in the Q. Therefore cross terms 
appear : 

total ~ .Normal + c t ,LIRC formal I LIRG 

+ 1-CcjJormallULIRG + CfilRG/ULIRc) (24) 



We plot in Fig. 16 the contributions of normal, LIRGs and 
ULIRGs. Note that the sum of the three contributions does not 
make the total power spectrum because the cross terms are not 
shown. 

Normal galaxies and LIRGs both dominate the power spectrum 
up to 550 fj.m. The contribution of LIRGs increases slightly 
and finally dominates from 850 yum to 2 mm. ULIRGs never 
clearly dominate the power spectrum at long wavelength, how- 
ever their relative contribution increases at long wavelengths, 
from 850 yum to 2 mm. 

Therefore, we do recover what is expected from previous works. 
Normal galaxies dominate at low redshift, LIRGs at z ~ 1 and 
ULIRGs contribute in the same way as the others at high redshift 
and thus at long wavelengths. 

6. Conclusion 

We presented a new model of the clustering of star-forming 
galaxies in the Cosmic Infrared Background anisotropies. We 
interfaced a parametric model of star-forming galaxies evolution 
with a halo distribution approach. The model is fully paramet- 
ric. Fixing the cosmology, it depends on the parameters of the 
model of galaxies and the HOD. We computed power spectra 
from 100 fim to 2 mm for IRAS, Spitzer/MIPS, Herschel/SPIRE 
and Planck/HFI spectral bands. We showed how power spectra 
can depend on the parameters and we concluded that the param- 
eters of the model of galaxies can hardly be constrained using 
C( only. Number counts and luminosity functions data are re- 
quired. Fixing them at the mean values found by Bethermin et al. 
( |201 1) , we explored the HOD parameters constraints and degen- 
eracies. The combination of Cc and counts/LF data do not break 
the degeneracies but constraints are slightly improved. Some of 
the parameters are strongly degenerate, especially M m ,„ and M sal 
with M mt = xM m i n with x ~ 3 where x is the direction of the de- 
generacy. 

We have shown that the 1-halo term can be detected at all wave- 
lengths and that galaxies at high redshift lie in smaller halos than 
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those at lower redshift. The level of the shot noise might not be 
reached with certain instruments such as Planck because of their 
angular resolution. However, this does not apply to the South 
Pole Telescope and to Herschel as they have a better angular res- 
olution. 

Using our fiducial model, we computed the halo mass and red- 
shift contribution to the power spectra. Higher redshift galax- 
ies contribute more at long wavelengths. Not surprisingly, the 
1- and 2- halo terms do not have the same mass dependence. We 
found that high mass halos contribute the most to the 1 halo term 
whereas the 2-halo term is dominated by intermediate mass ha- 
los which are most numerous. Our model strongly depends on 
the emissivity given by the evolution model of galaxies, and we 
compare the resulting Qs with those obtained using the emissiv- 
ities coming from Lagache et al. ( 2004). We have shown that the 
halo parameters strongly depend on the emissivities when data 
are fitted. 

In order to avoid the use of a model of evolution of galaxies, 
we have split the redshift range in four bins and computed the 
C( using the mean emissivity on these four redshift bins as in 
Amb lard et al.| ( p011| l. We considered these four parameters as 
free. They are not very well constrained, they cannot give any 
constraints on models of galaxies. 

We investigate the contribution of LIRGs and ULIRGs to the 
power spectra and its evolution with the wavelength. Our re- 
sults are in agreement with previous studies of normal galaxies, 
LIRGs and ULIRGs contribution to the CIB and to the lumi- 
nosity functions. Normal galaxies dominate the power spectrum 
at low redshift. As the redshift increases, LIRGs dominate Q. 
Meanwhile, the contribution of ULIRGs keeps increasing up to 
2 mm. 

The main unknown in CIB anisotropies power spectrum mea- 
surements are the redshift distributions of CIB galaxies. The 
coming results from Planck and Herschel will enable a great 
leap in the understanding of the clustering of star-forming galax- 
ies and its redshift evolution, measuring the cross power spec- 
tra between wavelengths. These new measurements will help to 
break some degeneracies and will allow to make more precise 
measurements of the star formation density at high redshift, and 
the characteristic mass of the dark matter halo at which the effi- 
ciency of the star formation is maximal. 
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Fig. 6. Counts at 160 fim, emissivities and clustering power spec- 
trum for three values of r^j z , the best fit and the best fit ±2<t. 
Pink crosses are data, Bether min et al.| ( |2010| ) for the counts and 
Lagache et al. ( 2007 1 for the power spectrum. Top panel : differ- 
ential number counts at 160 /mi. 2nd panel : mean emissivities at 
160 fim. 3rd panel : ratio of the modified emissivities compared 
to the best fit one at 160 /mi. 4th panel : power spectra of the 
clustering at 160 fim. Bottom panel : ratio of the modified power 
spectra compared to the best fit one at 160 fim. A small change 
of r Lir j z leads to a 20 % modification on the emissivities and on 
15-35% on the C t . 
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Fig. 10. Contribution of halo masses and redshift to the C( from 
100 |/m to 2 mm. The first column shows the redshift and mass 
contribution to the 1-halo term (£ = 2002) and the second 
column represents the same contributions to the 2-halo term 
({ = 100). The light grey corresponds to the highest contribu- 
tion to the Qs. The step of the color range is logarithmic and the 
scale is the same for both columns. The 1 -halo term is dominated 
by high masses at low z at all wavelength. At all wavelength, the 
2-halo term is dominated by a large range of masses at low z 
and by intermediate masses at higher z- The 1- and 2-halo term 
are sensitive to different mass regimes. The former to high mass 
halos and the latter to intermediate mass halos. 
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Fig. 11. Redshift distribution for several wavelengths for halo 
parameters different from the values of the fiducial model. The 
left panels are for logM m ,„ = 10.5 instead of log M m ;„ = 11.5 
and the right panels are for a — 1 instead of a sat = 1 .4. The red- 
shift distribution depends strongly on the chosen halo parameters 
although the same trend is observed : the low redshift dominate 
at short wavelengths, this contribution decreases with the wave- 
length and as the wavelength gets longer the contribution of high 
redshifts increases. 




Fig. 12. Linear bias for several values of k (in (Mpc//z) _1 ) and 
HOD bias. Light blue squares and diamonds bias values from 
resolved galaxies and pink crosses represent biases from un- 
resolved galaxies measurements (see Table [4|. The big blue 
squares represent the HOD biases coming from resolved galax- 
ies from Cooray et al. (2010). Light blue diamonds and squares 
show HOD and linear biases, respectively. 
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Fig. 13. Emissivities used by 



Viero et al. 



( 2009 1 coming from 



Lagache et al. ( 2004 1 (red continous line) and ours (black dot- 
ted line) coming from Bethermin et al. (2011 1 at several wave- 
lengths. 




Fig. 14. Power spectra ratios computed using the emissivities of 
Lagache et aT] ( |2004| i (C ejag ) and |Bethermin et al.| ( [20TT) (Q, MS ), 
for several wavelengths. 
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Fig. 15. Redshift contribution to the Q at several wavelengths using our emissivities on the left and using Lagache et al. (2004i 
emissivity on the right. In both cases, the contribution of high redshift increases with the wavelength but the evolution goes faster 
with Lagache et al. ( 2004| l emissivities. These redshift distributions are strongly varying with the halo parameters (see Fig. 10) 
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Fig. 16. Contribution to the Q of normal galaxies, LIRGs and ULIRGs at several wavelengths. When computing power spectra 
by splitting the contributions of normal galaxies, LIRGs and ULIRGs, cross terms appear. We do not not show them on that plot, 
therefore the sum of the three power spectra do not make the total. We recover that normal and LIRGs dominate at short wavelengths, 
thus at small redshifts. As the wavelength increases, and at the mean time the high-redshift contribution, the relative contribution of 
ULIRGs increases. 



